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ABSTRACT 


The  retrieval  of  an  aerosol  scattering  phase  function  using  a  multi-satellite 
technique  is  proposed.  A  total  of  33  phase  functions  were  derived  from  18 
smoke  cases  and  15  dust  cases.  Each  case  was  interrogated  using  four  to  nine 
satellite  passes  over  the  aerosol  region  in  a  two  to  four  hour  time  frame.  The 
radiance  values  for  the  red  and  near-infrared  (NIR)  channels  were  combined  with 
backscattering  angles  to  determine  the  shape  of  the  scattering  phase  function. 
The  radiance  values  were  input  into  the  Naval  Postgraduate  School  (NPS) 
aerosol  model  to  determine  optical  depths  and  sample  phase  functions.  A 
comparison  was  made  between  the  actual  phase  functions  retrieved  and  the 
NPS  model  phase  functions.  It  was  found  that  the  phase  functions  for  the  smoke 
cases  more  closely  matched  the  model  phase  functions  than  in  the  dust  cases. 
Some  conclusions  could  be  drawn  about  the  actual  aerosol  size  and  density 
distribution  based  on  how  well  it  matched  the  model  phase  function.  Further 
analysis  is  necessary  to  define  the  exact  size  and  number  of  the  aerosol 
particles.  Fully  understanding  the  aerosol  composition  is  crucial  in  determining 
its  effects  on  military  sensors  and  impacts  to  operations. 
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INTRODUCTION 


A.  MOTIVATION 

The  Air  Force  Weather  (AFW)  Strategic  Plan  and  Vision  outlines  eight 
“Transformation  Vectors”  that  will  radically  change  the  weather  core  processes  of 
collection,  analysis,  prediction,  tailoring,  and  dissemination.  Transformation 
vector  #1  is  to  integrate  environmental  information  collection  into  the  theater 
Intelligence,  Surveillance  and  Reconnaissance  (ISR)  collection  plan.  Issue  #1  of 
transformation  vector  #1  is  to  improve  the  ability  to  collect  globally  target-quality 
environmental  information  in  hostile  and  data-denied  areas  to  build  an  accurate 
real-time  picture  of  the  environment  for  battlespace  awareness  and  decision 
superiority  (Air  Force  Weather,  2004).  Obtaining  a  complete  picture  of 
atmospheric  aerosol  properties  is  a  crucial  step  in  obtaining  full  battlespace 
awareness. 

Satellite  data  from  several  polar  orbiting  satellites  was  used  in  this  thesis 
to  determine  the  scattering  phase  function  of  an  aerosol  region.  Data  from  the 
Moderate  Resolution  Imaging  Spectroradiometer  (MODIS)  instrument  onboard 
the  Earth  Observing  System  (EOS)  satellites,  the  Advanced  Very  High 
Resolution  Radiometer  (AVFIRR)  onboard  the  National  Oceanographic  and 
Atmospheric  Administration  (NOAA)  satellites,  and  the  Sea-viewing  Wide  Field- 
of-view  Sensor  (SeaWiFS)  onboard  the  National  Aeronautics  and  Space 
Administration  (NASA)  satellite  was  used  to  determine  the  scattering  phase 
function  for  various  aerosol  areas. 

Currently  NASA's  Multi-angle  Imaging  SpectroRadiometer  (MISR)  (Figure 
1)  is  performing  a  similar  technique  for  aerosol  and  terrain  interrogation.  The 
MISR  views  the  earth  in  a  narrow  swath  from  nine  different  view  angles.  Most 
satellite  instruments  look  only  straight  down,  or  toward  the  edge  of  the  planet.  To 
fully  understand  Earth's  climate,  and  to  determine  how  it  may  be  changing,  we 
need  to  know  the  amount  of  sunlight  that  is  scattered  in  different  directions  under 
natural  conditions.  MISR  is  an  instrument  designed  to  address  this  need;  it  views 
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the  Earth  with  sensors  pointed  at  nine  different  angles.  One  sensor  points  toward 
nadir,  and  the  others  provide  forward  and  aftward  view  angles,  at  the  Earth's 
surface,  of  26.1°,  45.6°,  60.0°,  and  70.5°.  As  the  instrument  flies  overhead,  each 
region  of  the  Earth's  surface  is  successively  imaged  by  all  nine  sensors  in  each 
of  four  wavelengths  (blue,  green,  red,  and  near-infrared).  With  its  nine  different 
look  angles,  MISR  is  uniquely  capable  of  observing  cloud  and  aerosol  plume 
structures  in  the  atmosphere  three-dimensionally,  as  well  as  measuring  the 
relative  heights  and  types  of  clouds.  Its  stereoscopic  measurement  abilities 
enable  MISR  to  observe  the  amount  of  sunlight  that  is  reflected  at  various 
angles.  MISR  will  also  help  scientists  trace  smoke  and  aerosol  plumes  back  to 
their  sources  (NASA  Jet  Propulsion  Laboratory  2005). 


Figure  1 .  Multi-angle  Imaging  SpectroRadiometer  (MISR) 

B.  OBJECTIVES 

Improving  visibility  estimates  in  areas  of  limited  data  has  been  and  will 
continue  to  be  an  area  of  interest  for  the  military  war-fighter.  Many  tools  are 
available  to  assist  in  the  visibility  estimates.  However,  most  are  based  on  a 
single  sized  spherical  aerosol  model.  Unfortunately  aerosols  are  often  mixed  in 
composition,  size  and  shape.  This  leads  to  a  very  complex  scattering  regime  of 
the  aerosols,  which  is  not  handled  well  by  any  single  aerosol  model.  The 
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purpose  of  this  thesis  is  to  determine  the  scattering  phase  function  of  an  aerosol 
of  unknown  size,  shape,  and  composition.  This  in  turn  would  lead  to  better 
optical  depth  interpretation,  which  would  lead  to  better  visibility  estimation. 

The  thesis  will  start  with  a  description  of  the  instruments  used  for  data 
collection.  That  is  followed  by  a  brief  explanation  of  scattering  theory.  A  large 
section  is  devoted  to  the  procedures  used  to  determine  the  scattering  phase 
function.  The  results  section  describes  the  findings  of  this  study.  The  last 
section  contains  the  conclusions  and  some  recommendations  for  further  study. 


3 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


4 


II.  BACKGROUND 


A.  SATELLITE  INSTRUMENTS 

Several  different  satellite  instruments  were  used  for  the  collection  of 
radiance  values  in  the  red  and  near-infrared  (NIR)  wavelengths 

1.  Advanced  Very  High  Resolution  Radiometer  (AVHRR) 

The  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  provides  as 
few  as  four  to  as  many  as  six  bands  of  multispectral  data  from  the  NOAA  polar- 
orbiting  satellite  series.  There  has  been  fairly  continuous  global  coverage  since 
June  1979,  with  morning  and  afternoon  acquisitions  available  (Table  1).  The 
resolution  is  1 .1  kilometer  at  nadir. 


Satellite  # 

Dates 

Pass 

#  of  bands 

NOAA-6 

6/79-11/86 

am 

4 

NOAA-7 

8/81  -  6/86 

pm 

5 

NOAA-8 

5/83-10/85 

am 

4 

NOAA-9 

2/85-11/88 

pm 

5 

NOAA- 10 

11/86-9/91 

am 

4 

NOAA-1 1 

11/88-9/94 

pm 

5 

NOAA- 12 

5/91  -12/94 

am 

5 

NOAA-1 3 

failed 

pm 

n/a 

NOAA-1 4 

1 2/94  -  present 

pm 

5 

NOAA-1 5 

5/98  -  present 

am 

5 

NOAA-1 6 

9/00  -  present 

pm 

6 

NOAA-1 7 

1 0/02  -  present 

am 

6 

Table  1 .  NOAA  satellite  properties  and  history 
(from  US  Geological  Survey,  2004) 

In  Table  2  the  NOAA  bands  are  listed.  Bands  one  and  two  were  used  to 
collect  the  red  and  NIR  data  respectively.  The  spectral  response  curves  for 
Bands  one,  two  and  three  are  shown  in  Figures  2,  3  and  4.  Their  significance  is 
explained  in  further  detail  in  the  Procedures  section  of  this  thesis. 
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Band 

# 

Satellites: 

NOAA- 

Satellites: 

NOAA- 

Satellites: 
NOAA-1 5,1 6,1 7 

IFOV 

6,8,10 

7,9,11,12,14 

1 

0.58-0.68 

0.58-0.68 

0.58-0.68 

1.39 

2 

0.725-1.10 

0.725-1.10 

0.725-1.0 

1.41 

3 

3.55-3.93 

3.55-3.93 

1.51 

3A 

1.58-1.64 

1.51 

3B 

3.55-3.93 

1.51 

4 

10.50-11.50 

10.3-11.3 

10.3-11.3 

1.41 

5 

band  4  repeated 

11.5-12.5 

11.5-12.5 

1.30 

(micrometers) 

(micrometers) 

(micrometers) 

(milliradians) 

Table  2.  NOAA  AVHRR  Bands 

(from  US  Geological  Survey,  2004) 
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Figure  2.  Spectral  Response  Curve  for  NOAA-17  AVHRR/3  Channel  1 

(from  NOAA,  2002) 
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Figure  3.  Spectral  Response  Curve  for  NOAA-1 7  AVHRR/3  Channel  2 

(from  NOAA,  2002) 
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Figure  4.  Spectral  Response  Curve  for  NOAA-1 7  AVFIRR/3  Channel  3A/3B 

(from  NOAA,  2002) 


2.  Moderate  Resolution  Imaging  Spectroradiometer  (MODIS) 
MODIS  (or  Moderate  Resolution  Imaging  Spectroradiometer)  is  a  key 
instrument  aboard  the  Terra  (EOS  AM)  and  Aqua  (EOS  PM)  satellites.  Terra’s 
orbit  around  the  Earth  is  timed  so  that  it  passes  from  north  to  south  across  the 
equator  in  the  morning  (10:30  a.m.  LT),  while  Aqua  passes  south  to  north  over 
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the  equator  in  the  afternoon  (1 :30  p.m.  LT).  Terra  MODIS  and  Aqua  MODIS  are 
viewing  the  entire  earth’s  surface  every  1  to  2  days,  acquiring  data  in  36  spectral 
bands,  or  groups  of  wavelengths  (Tables  3  and  4).  Bands  one  and  two  were 
used  to  collect  the  red  and  NIR  data  respectfully. 


Prim  ary  Use 

Band 

Bandw  i d t h 1 

Spectral 

Radiance 

Required 

SNRa 

Land  /  Cloud  / Aerosols 
Boundaries 

1 

620  -  670 

21 .8 

128 

2 

841  -  876 

24.7 

201 

Land  /  Cloud/ Aerosols 
Properties 

3 

459  -  479 

35.3 

243 

4 

545  -  565 

29.0 

228 

5 

1230  -  1250 

5.4 

74 

6 

1628  -  1652 

7.3 

275 

7 

2105  -  2155 

1.0 

110 

Ocean 

Color/ 

Phytoplankton/ 
Biogeochem  istry 

8 

405  -  420 

44.9 

880 

9 

438  -  448 

41 .9 

838 

10 

483  -  493 

32.1 

802 

1 1 

526  -  536 

27.9 

754 

12 

546  -  556 

21.0 

750 

13 

662  -  672 

9.5 

910 

14 

673  -  683 

8.7 

1087 

15 

743  -  753 

10.2 

586 

16 

862  -  877 

6.2 

516 

Atm  ospheric 

Water  Vapor 

17 

890  -  920 

10.0 

167 

18 

931  -  941 

3.6 

57 

19 

915  -  965 

15.0 

250 

Table  3.  MODIS  Bands  1-19 
(from  Conboy,  2002) 
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Prim  ary  Use 

Band 

Bandw  idth 

Spectral 

Radiance- 

Required 
N  E[  delta]  T(  K)1 

Surface/  Cloud 
Tern  perature 

20 

3.660  -  3.840 

0.45(300K) 

0.05 

21 

3.929  -  3.989 

2.38(335K) 

2.00 

22 

3.929  -  3.989 

0.67(300K) 

0.07 

23 

4.020  -  4.080 

0. 79(300K) 

0.07 

Atm  ospheric 
Tern  perature 

24 

4.433  -  4.498 

0. 1  7(250K) 

0.25 

25 

4.482  -  4.549 

0.59(275K) 

0.25 

Cirrus  Clouds 
Water  Vapor 

26 

1.360  -  1.390 

6.00 

1  50(SNR) 

27 

6.535  -  6.895 

1 . 1  6(240K) 

0.25 

28 

7.175  -  7.475 

2. 1  8(250K) 

0.25 

Cloud 

Properties 

29 

8.400  -  8.700 

9.58(300K) 

0.05 

Ozone 

30 

9.580  -  9.880 

3. 69(250K) 

0.25 

Surface/  Cloud 
Tern  perature 

31 

10.780  -  11.280 

9.55(300K) 

0.05 

32 

1  1 .770  -  12.270 

8. 94(300K) 

0.05 

Cloud  Top 
Altitude 

33 

13.185  -  13.485 

4.52(260K) 

0.25 

34 

13.485  -  13.785 

3. 76(250K) 

0.25 

35 

13.785  -  14.085 

3.11  (240K) 

0.25 

36 

14.085  -  14.385 

2. 08(220K) 

0.35 

1  Bands  1  to  1  9  are  in  nm;  Bands  20  to  36  are  in  pm 

2  Spectral  Radiance  values  are  (W/m2  -pm-sr) 

3SNR=  Signal-to-noise  ratio 

4  NE(delta)T  =  Noise-equivalent  temperature  difference 

Note:  Performance  goal  is  30-40%  better  than  required 

Table  4.  MODIS  Bands  20-36 
(from  Conboy,  2002) 


3.  Sea-viewing  Wide  Field-of-view  Sensor  (SeaWiFS) 

The  purpose  of  the  Sea-viewing  Wide  Field-of-view  Sensor  (SeaWiFS) 
project  is  to  provide  quantitative  data  on  global  ocean  bio-optical  properties  to  the 
Earth  science  community.  The  NASA  SeaWiFS  instrument  is  composed  of  an 
optical  scanner  and  an  electronics  module  (Figure  5).  The  eight  bands  (Table  5), 
most  of  which  are  visible  wavelengths,  make  it  very  useful  for  the  interrogation  of 
aerosols.  Another  key  feature  of  the  SeaWiFS  sensor  is  its  scanner  tilt 
mechanism  that  enables  the  instrument  to  be  oriented  in  the  along-track  direction 
from  +20  to  -20  degrees  to  avoid  sun  glint  from  the  sea  surface.  Tilting  the  entire 
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scanner,  rather  than  only  a  section  of  the  optical  train,  assures  that  the  SeaWiFS 
calibration,  polarization,  and  angular  scanning  characteristics  will  be  identical  for 
all  tilt  positions  and,  thereby,  simplifies  the  ground  processing  of  in-flight  data. 


Instrument  Bands 

Band 

Wavelength 

1 

402-422  nm 

2 

433-453  nm 

3 

480-500  nm 

4 

500-520  nm 

5 

545-565  nm 

6 

660-680  nm 

7 

745-785  nm 

8 

845-885  nm 

Mission  Characteristics 

Orbit  Type 

Sun  Synchronous  at  705  km 

Equator  Crossing 

Noon  +20  min,  descending 

Orbital  Period 

99  minutes 

Swath  W  idth 

2,801  km  LAC/HRPT  (58.3  degrees) 

Swath  W  idth 

1,502  km  GAC  (45  degrees) 

Spatial  Resolution 

1.1  km  LAC,  4.5  km  GAC 

Real-Time  Data  Rate 

665  kbps 

Revisit  Time 

1  day 

Digitization 

1 0  bits 

Table  5.  SeaWiFS  Bands  and  Mission  Characteristics 
(from  Feldman,  1997) 
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Figure  5.  NASA  Satellite  and  SeaWiFS  Sensor 
(from  Feldman,  1997) 


B.  AERONET 

(AErosol  RObotic  NETwork)  is  an  optical  ground  based  aerosol  monitoring 
network  and  data  archive  supported  by  NASA's  Earth  Observing  System  and 
expanded  by  federation  with  many  non-NASA  institutions.  The  network 
hardware  consists  of  identical  automatic  sun-sky  scanning  spectral  radiometers 
owned  by  national  agencies  and  universities.  Data  from  this  collaboration 
provides  globally  distributed  near  real  time  observations  of  aerosol  spectral 
optical  depths,  aerosol  size  distributions,  and  precipitable  water  in  diverse 
aerosol  regimes.  The  data  undergo  preliminary  processing  (real  time  data), 
reprocessing  (final  calibration  ~6  mo.  after  data  collection),  quality  assurance, 
archiving  and  distribution  from  NASA's  Goddard  Space  Flight  Center  master 
archive  and  several  identical  data  bases  maintained  globally.  The  data  provide 
algorithm  validation  of  satellite  aerosol  retrievals,  as  well  as  characterization  of 
aerosol  properties  that  are  unavailable  from  satellite  sensors. 

The  primary  instrument  of  the  AERONET  is  the  CIMEL  Electronique  31 8A 
spectral  radiometer  (Figure  6)  is  a  solar  powered  weather  hardy  robotically 
pointed  sun  and  sky  spectral  radiometer.  The  radiometer  makes  two  basic 
measurements,  either  direct  sun  or  sky,  both  within  several  programmed 
sequences.  The  direct  sun  measurements  are  made  in  eight  spectral  bands  at 
wavelengths  of  340,  380,  440,  500,  670,  870,  940  and  1020  nm.  A 
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preprogrammed  sequence  of  measurements  is  taken  by  these  instruments. 
Optical  thickness  is  calculated  from  spectral  extinction  of  direct  beam  radiation  at 
each  wavelength  based  on  the  Beer-Bouguer  Law.  Attenuation  due  to  Rayleigh 
scatter,  absorption  by  ozone  (from  interpolated  ozone  climatology  atlas),  and 
gaseous  pollutants  is  estimated  and  removed  to  isolate  the  aerosol  optical 
thickness  (AOT).  Sky  radiance  measurements  are  inverted  with  the  Dubovik  and 
Nakajima  inversions  to  provide  aerosol  properties  of  size  distribution  and  phase 
function  over  the  particle  size  range  of  0.1  to  5  pm  (Holbren,  2004).  Although  the 
data  from  the  AERONET  is  highly  accurate  it  has  limited  coverage.  It  was  often 
difficult  to  locate  a  sensor  in  close  proximity  to  the  over  the  water  aerosol  regions 
being  interrogated  .Displayed  in  Figure  7  is  a  sample  of  the  AERONET  optical 
thickness  output.  All  of  the  AERONET  data  products  used  can  be  seen  in 
Appendix  C.  The  graphics  were  generated  at  the  NASA,  Goddard  Space  Flight 
Center  web  site  (http://aeronet.gsfc.nasa.gov/data_menu.html). 


Figure  6.  CIMEL  Spectral  Radiometer 
(from  Holbren,  2004) 
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Figure  7.  AERONET  Optical  Thickness  output  for  Moldova  near  the  Black  Sea 

(from  NASA  GSFC,  2005a) 

C.  SCATTERING 

Radiation  scatted  from  a  particle  is  a  function  of  particle  size,  particle 
shape,  particle  index  of  refraction,  wavelength  of  radiation,  and  viewing 
geometry.  For  a  spherical  scatterer,  the  scattered  radiation  is  a  function  of  only 
the  viewing  angle,  index  of  refraction,  and  the  size  parameter  defined  as: 


__  2 Jtr_ 

~  a  (i) 

where  r  is  the  radius  of  the  sphere  (from  Kidder  and  Vonder  Flaar,  1995).  The 
size  parameter  can  be  used  to  divide  scattering  into  three  regimes:  Geometric, 
Rayleigh,  and  Mie  (Figure  8).  The  latter  being  the  regime  of  interest  for  aerosols 
in  the  visible  and  near-infrared  wavelengths. 
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Figure  8.  Scattering  Regimes 

(from  Wallace  and  Hobbs,  1977) 

1.  Geometric 

For  X  greater  than  about  50,  the  sphere  is  large  compared  to  the 
wavelength  of  radiation.  In  this  regime,  rays  which  are  reflected  and  refracted  at 
the  surface  of  a  scatterer  can  be  traced.  Sun  glint  reflected  off  of  the  surface  of 
the  water  is  an  example  of  geometric  optics.  The  water  surface  is  somewhat 
rough  therefore  the  solar  reflection  is  blurred  and  larger  than  the  sun. 

2.  Rayleigh 

For  X  less  than  about  0.1,  the  particle  is  small  compared  to  the 
wavelength  of  radiation.  In  this  regime,  atmospheric  molecules  act  as  Rayleigh 
scatterers  for  visible  and  ultraviolet  radiation.  Rayleigh  scattering  is  insensitive  to 
particle  shape.  Rayleigh  particles  scatter  equally  well  in  the  forward  as  well  as 

backward  direction  (Figure  9).  Scattering  is  negligible  for  X  less  than  10'3. 

3.  Mie 

For  X  in  the  range  of  0.1-50,  the  particle  size  is  comparable  to  the 

wavelength  of  radiation.  It  is  in  the  Mie  regime  that  aerosols  interact  with  the 
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visible  and  near-infrared  wavelengths.  The  scattering  properties  become  very 
complicated  in  the  Mie  regime.  The  scattering  phase  function  is  dependent  upon 
the  angle  of  incident  radiation  to  the  aerosol.  It  can  be  seen  in  Figure  9  that 
there  is  a  very  strong  forward  scattering  component  that  increases  logarithmically 

as  the  ^-value  increases.  It  is  also  noted  that  the  backscattered  radiation  varies 
dramatically  with  the  backscattering  angle. 


Figure  9.  Polar  plots  of  the  scattering  phase  function  for  several  size  parameters 

(from  Kidder  and  Vonder  Flaar,  1995) 

D.  ASSUMPTIONS 

Some  necessary  assumptions  had  to  be  made  in  working  with  the  satellite 
data.  The  total  radiance  reaching  the  satellite  sensor  can  be  broken  down  into 
the  following  expression  (from  Durkee  et  al. ,  1986): 

Lt  =  LoT  o  +  Lr  +  La  (2) 

where  Lt  is  the  total  radiance  reaching  the  sensor,  L0  is  the  radiance  leaving  the 
surface,  T0  is  the  transmittance  from  the  surface  to  the  sensor,  Lr  is  the  radiance 
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from  Rayleigh  scattering  by  atmospheric  molecules,  and  La  is  the  radiance 
backscattered  by  the  aerosol.  An  illustration  of  the  radiance  paths  reaching  the 
satellite  can  be  seen  in  Figure  10. 


^  “  Lr  +  *"A  +  0-G  +  ^)rd 


Lg  =  Surface  Radiance 
=  Glint  Radiance 
Lr  =  Rayleigh  Radiance 
La  =  Aerosol  Radiance 

Uo  =  l-s  +  Ug 


Molecular 

Rayleigh 

Scattering 


Aerosol 
Mie 

Scattering 

Surface  Geometric 
Scattering 


t  0)0  rr  /  \^a 

La  =  —  Eop(y/)  — 

An  fi 


over  ocean  surface 

Figure  10. 


Breakdown  of  total  radiance 
(after  Durkee,  2004) 


It  is  assumed  that  the  molecular  Rayleigh  scattering  remains  constant 
during  the  sample  period.  It  can  be  seen  in  Figure  9  that  the  Rayleigh  scattering 
is  more  nearly  constant  in  all  directions  then  for  Mie  scattering.  It  is  also 
assumed  that  the  surface  radiance  and  transmittance  do  not  change  over  the 
sample  period.  With  these  assumptions  that  leaves  the  radiance  from  the 
aerosol  as  the  dominant  contributor  to  the  changes  in  the  total  radiance. 
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The  procedures  section  that  follows  describes  how  the  total  radiance 
values  were  determined  from  the  satellite  data.  The  total  radiance  values  vary 
form  one  satellite  to  the  next  because  the  Mie  scattering  from  the  aerosol  is 
changing  as  the  satellite  and  solar  angles  are  changing.  The  procedures  section 
describes  the  complete  process  for  determining  the  backscattering  angle  and 
combining  it  with  the  radiance  values  to  ultimately  produce  a  scattering  phase 
function  for  the  aerosol. 
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III.  PROCEDURES 


A.  CASE  SELECTION 

Cases  were  selected  based  on  the  satellites’  ability  to  visibly  see  aerosols 
and  track  their  movement  over  water.  Higher  latitudes  are  better,  since  this 
increased  the  chance  of  getting  overlapping  passes  on  the  subsequent  orbit  from 
a  single  satellite.  The  selection  of  a  case  required  there  to  be  a  cloud  free 
environment.  Clouds  would  contaminate  the  radiance  values  backscattered  from 
the  aerosols.  Sun  glint  also  caused  a  problem  with  the  radiance  that  reached  the 
satellite  sensor.  It  was  therefore  necessary  to  ensure  that  no  sun  glint  was  in  the 
area  of  the  aerosols  being  interrogated.  Once  a  possible  case  was  identified  the 
raw  satellite  data  was  ordered  from  archives  and  then  loaded  into  Terascan  for 
viewing  and  processing. 

B.  TERASCAN 

Terascan  was  used  to  display  the  raw  satellite  data  from  NOAA,  EOS,  and 
SeaWiFS.  The  Terascan  software  allowed  for  the  retrieval  of  all  pertinent  data  at 
any  specific  pixel.  The  values  of  interest  for  these  cases  were  radiance  values  in 
the  red-visible  and  near-IR  wavelengths,  solar  zenith  angle,  solar  azimuth  angle, 
sensor  zenith  angle,  and  sensor  azimuth  angle.  The  radiance  values  for  the 
aerosol  were  determined  by  sampling  a  100-pixel  area  over  the  aerosol  region. 
Standard  statistical  methods  from  Terascan  were  used  to  determine  the  mean 
and  standard  deviation  of  the  radiance  values  in  the  sample  area.  This  allowed 
for  greater  confidence  in  using  the  radiance  value  to  determine  the  scattering 
phase  function.  The  values  at  the  center  of  the  100-pixel  box  were  used  for  solar 
zenith,  solar  azimuth,  sensor  zenith,  and  sensor  azimuth.  It  was  necessary  to 
obtain  these  radiance  and  angle  values  from  a  series  of  satellite  images,  each  of 
which  viewed  the  aerosol  region  of  interest  from  different  view  angles. 

C.  DATA  RETRIEVAL  AND  INTEROGATION 

1.  Radiance  Values 

For  the  red  wavelength:  Channel  1  (0.620-0.670 nm)  was  used  for 
MODIS,  Channel  1  (0.580-0.680 jum)  was  used  for  AVHRR,  and  Channel  6 
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tj  Survey  Report 


modi s  chOl 


Minimum  ; 

Maximum  t 
Mean  : 

Standard  Deviation  : 
Number  of  Samples  ; 
Humber  of  Good  Samples 


(0.660-0.680 pm)  was  used  for  SeaWiFS.  For  the  near-infrared  wavelength: 
Channel  2  (0.841  -0.876 pm)  was  used  for  MODIS,  Channel  2  (0.725-1.0 pm) 
was  used  for  AVFIRR,  and  Channel  8  (0.845-0.885 pm )  was  used  for  SeaWiFS. 

To  determine  the  radiance  values  a  100-pixel  sample  was  taken  within  the 
aerosol-laden  regions  of  the  scene.  The  radiance  value  for  each  channel  was 
taken  as  the  mean  value  of  that  1 00-pixel  sample  (Figure  1 1 ). 


Figure  1 1 .  Statistical  Sample  of  radiances  for  MODIS  Channel  1 . 


Various  conversion  factors  were  used  to  ensure  all  radiance  outputs  were 
in  the  same  units  (watts/m2/ pm).  The  MODIS  data  required  no  conversion  while 


alue 

I  nodis_ch01 


-  V  Tine  03/09/09  10;45;05 _ 

J  Var  |  nodis_ch01  r~ 
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the  AVHRR  and  SeaWiFS  data  required  some  conversions.  The  AVHRR  output 
was  in  albedo  which  required  the  following  conversion: 

(Albedo)*E0/^  =  (watts/m2/  fim),  (3) 

where  E0  is  the  solar  constant  at  the  centered  effective  wavelength.  The  values 
of  E0  for  AVHRR  channel  1  and  2  were  1644  watts/m2/ fim  and  1034 
watts/m2/ fim  respectively.  The  output  for  SeaWiFS  was  in  radiance  and  only 
required  a  correction  factor  of  10  to  be  in  the  standard  units: 

(SeaWiFS  mW/cm2/ fim  )*1 0  =  (watts/m2/  fim )  (4) 

2.  Linear  Extrapolation 

A  linear  extrapolation  was  necessary  to  ensure  that  all  the  radiance  values 
were  being  compared  at  the  same  wavelength.  The  wavelengths  for  the  AVHRR 
on  NOAA  14  were  used  in  order  to  match  the  NPS  aerosol  model  described 
below.  For  the  red  wavelength:  MODIS  Channel  1  has  a  centered  effective 
wavelength  at  0.645  fim ;  AVHRR  Channel  1  has  a  centered  effective  wavelength 
at  0.6357//m;  and  SeaWiFS  Channel  6  has  a  centered  effective  wavelength  at 
0.670  fim .  For  the  near-infrared  wavelength:  MODIS  Channel  2  has  a  centered 
effective  wavelength  at  0.8585 fim  ;  AVHRR  Channel  2  has  a  centered  effective 
wavelength  at  0.85397 fim  ]  and  SeaWiFS  Channel  8  has  a  centered  effective 
wavelength  at  0.865  fim  (NOAA,  2002).  A  linear  extrapolation  was  accomplished 
for  MODIS  channels  1  and  2  by  finding  the  slope  of  the  radiance  versus 
wavelength  line  between  channels  1  and  2,  then  applying  a  linear  extrapolation 
to  the  AVHRR  center  effective  wavelength  of  0.6357 fim  for  the  red  wavelength 
and  0.85397  fim  for  the  near-infrared  wavelength. 

The  following  formula  is  an  example  of  the  linear  extrapolation  of  the 
MODIS  channel  1  radiance  value  to  the  radiance  value  at  the  AVHRR  Channel  1 
center  effective  wavelength: 


Lac  = 


L2  —  L\ 
(A2  —  A1) 


(Ai  —  A  ac )  +  Li 


(5) 
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where  Li  and  L2  are  the  radiance  values  for  MODIS  channels  1  and  2;  k\  and 
/liare  the  wavelengths  for  MODIS  channels  1  and  2;  A,ac  is  the  center  effective 
wavelength  for  AVHRR  channel  1;  and  LAc  is  the  linear  extrapolated  radiance 
value  for  the  red  wavelength.  A  graphic  illustration  of  the  linear  extrapolation  can 
be  seen  in  Figure  12. 
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Figure  12.  Linear  Extrapolation  of  radiance  values  to  AVHRR  red  and  NIR  center 

effective  wavelengths 


The  same  procedures  and  similar  formulas  were  used  to  extrapolate  the 
radiance  values  for  MODIS  channel  2,  SeaWiFS  channels  6  and  8  to  the 
radiance  values  of  the  center  effective  wavelengths  of  AVHRR  channels  1  and  2. 
The  net  result  is  that  all  radiance  values  were  extrapolated  to  radiance  values 
centered  at  0.6357 jum  for  the  red  wavelength  and  0.85397 jum  for  the  near- 
infrared  wavelength. 

3.  Phase  Scattering  Angle 

The  solar  zenith  (O'),  solar  azimuth  {</>’),  sensor  zenith  (6),  and  sensor 
azimuth  {</>)  angles  obtained  from  Terrascan  were  used  to  determine  the 
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backscattering  angle  (y/).  The  values  for  the  angles  were  taken  from  the  center 
of  the  100-pixel  sample  box  and  used  in  the  following  equations  for  the  MODIS 
data. 


cos  if/'-  cos0cos0'+sin0sin#'cos(^-^')  (6) 

^  =  180-^'  (7) 

The  backscattering  angle  {y/)  was  obtained  directly  from  equation  seven 
for  the  AVHRR  and  SeaWiFS  data,  since  the  cosy/'  was  one  of  the  outputs  in 
terascan. 

4.  Plotting  the  Phase  Scattering  Function 

The  backscattering  angles  were  used  with  the  radiance  values  in  plotting 
the  scattering  phase  function.  The  following  formula  was  used  to  convert  the 
radiance  values  into  a  phase  function  as  a  function  of  the  backscattering  angle 
(from  Durkee  et  al.,  1986). 


La 


CQo 

4  n 


Eop{y/) 


Sa 


(8) 


and, 


La/dAn 

EoSaCOo  ’  (9) 

where  La  is  the  aerosol  radiance  at  the  center  effective  wavelength,  n  is  the 
solar  zenith  angle,  E0  is  the  aerosol  solar  constant  at  the  center  effective 

wavelength,  Sa  is  the  aerosol  optical  depth,  and  coQ  is  the  single  scattering 
albedo.  All  of  the  values  at  this  point  are  known  with  the  exception  of  the  optical 
depth. 

The  aerosol  optical  depth  was  determined  using  the  NPS  aerosol  model 
(Figure  13).  The  optical  depths  from  the  NPS  model  were  compared  against  any 
AERONET  observations  in  close  proximity.  This  was  to  ensure  the  optical  depth 
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values  from  the  model  were  a  reasonable  comparison  to  the  values  observed 
throughout  the  day  by  the  AERONET  sensor.  The  AERONET  optical  depth  data 
can  be  seen  in  Appendix  C. 
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Figure  13.  NPS  Aerosol  Model 
(from  Durkee,  2004) 


The  NPS  aerosol  model  starts  with  an  assumed  size  and  number 
distribution  of  the  aerosol  (Upper  right  of  Figure  13).  Based  on  that  size  and 
number  distribution  the  ratio  of  the  red  to  NIR  aerosol  radiances  is  determined. 
This  is  referred  to  as  the  Si2  value.  Based  on  the  Si2  value,  the  satellite  angles 
and  the  solar  angles,  an  aerosol  model  number  is  picked.  Each  aerosol  model 
number  has  a  predetermined  phase  function.  The  shape  of  the  phase  function  is 
derived  from  the  assumed  aerosol  size  and  number  distribution. 
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In  this  thesis  the  calculated  aerosol  scattering  phase  function  was 
retrieved  by  plotting  the  values  of  Equation  (9)  against  its  corresponding 
backscattering  angle  for  the  various  satellite  views  of  the  aerosol  region.  A 
comparison  of  the  calculated  and  model  phase  functions  is  done  in  the  results 
section.  The  results  section  also  contains  the  NPS  model  outputs  for  the  18 
smoke  aerosol  and  15  dust  aerosol  cases. 
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IV.  RESULTS 


A.  SMOKE 

A  total  of  18  smoke  cases  were  observed.  Each  case  had  between  four 
and  seven  satellite  views  over  the  aerosol  region  within  a  two  to  four  hour  time 
period.  Each  specific  case  can  be  seen  in  detail  in  Appendix  A.  The  smoke 
cases  were  run  on  the  NPS  aerosol  model  and  the  S12  ratio  computed.  It  was 
observed  that  the  SeaWiFS  sensor  often  produced  the  lowest  S12  value  (Figure 
14).  The  tight  grouping  of  the  S12  values  gives  confidence  in  the  multiple  satellite 
technique.  The  tight  grouping  indicates  that  the  satellites  are  viewing  similar 
ratios  of  the  red/NIR  wavelengths.  The  tightest  groupings  occurred  in  cases  1-3 
and  15-17,  when  the  smoke  optical  depths  were  highest  (Figures  15  and  16). 
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Figure  14.  S12  Ratio  for  Smoke  Cases 
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Figure  15.  NPS  Aerosol  Model  Optical  Depth  (AOD63)  for  the  Smoke  Cases 

The  NPS  aerosol  model  optical  depth  values  for  the  red  and  NIR 
wavelengths  are  indicated  in  the  AOD63  and  AOD86  Figures  15  and  16 
respectfully.  The  variability  in  the  optical  depth  seen  in  Figure  15  is  largest  when 
the  optical  depth  is  higher.  It  is  the  most  evident  in  the  NIR  wavelength  (Figure 
16)  where  the  lower  optical  depths  vary  only  slightly.  The  model  optical  depths 
are  within  the  range  of  optical  depth  values  observed,  in  a  24  hour  period,  by 
AERONET  in-situ  sensors. 
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NPS  Aerosol  Model  (AOD86)--Smoke 
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Figure  16.  NPS  Aerosol  Model  Optical  Depth  (AOD86)  for  the  Smoke  Cases 
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The  NPS  aerosol  model  number  is  chosen  based  on  what  the  model 
interprets  the  size  of  the  aerosol  particles  to  be.  The  smaller  particles,  like 
smoke,  have  a  higher  S12  value.  A  higher  S12  value  results  in  a  higher  model 
number.  The  NPS  aerosol  model  most  frequently  picks  model  the  highest  model 
number,  seven,  for  smoke  cases.  This  is  expected  since  model  seven  has  the 
smallest  particle  size  distribution.  The  model  numbers  shown  in  Figure  17 
indicate  that  the  model  number  from  the  SeaWiFS  sensor  is  frequently  lower 
than  that  of  the  other  sensors.  This  could  be  an  indication  of  a  calibration  issue 
for  the  SeaWiFS  sensor. 
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Figure  1 7.  NPS  Aerosol  Model  Phase  Function  Pick  for  the  1 8  Smoke  Cases 
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Figure  18  illustrates  the  retrieved  phase  function  values  from  all  satellites 
for  the  three  smoke  cases  over  Hudson  Bay.  It  is  apparent  that  the  phase 
function  is  decreasing  as  the  backscatter  angle  goes  from  100° -140°  and 
increasing  for  backscatter  angles  above  140°.  As  expected  the  calculated 
values  for  the  phase  function  are  higher  than  those  of  the  model.  The  calculated 
values  used  total  radiance  instead  of  the  aerosol  radiance.  Since  the  NPS  model 
removes  Rayleigh  radiance  and  makes  an  ozone  correction,  its  values  are 
smaller.  This  difference  is  not  important  because  the  primary  interest  in  this 
study  is  the  comparison  of  the  shape  of  the  calculated  vs.  model  phase  functions. 
The  phase  functions  vary  somewhat  in  shape  but  trend  toward  increasing  values 
after  140°.  The  comparative  shape  of  the  plotted  phase  function  to  that  of  the 
model  leads  to  the  conclusion  that  the  actual  size  distribution  of  the  aerosol 
particles  is  close  to  that  indicated  by  the  model.  However,  there  are  definitely 
some  differences  that  indicate  the  actual  size  and  number  of  aerosol  particles  are 
different  from  that  of  the  model  size  and  number  distribution. 
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Figure  18.  Phase  Function  for  Smoke  Cases  1-3,  over  Hudson  Bay,  with  NPS 

Aerosol  Model  #7  Phase  Function 
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The  retrieved  phase  functions  for  four  smoke  cases  over  the  Black  Sea 
are  combined  and  displayed  in  Figure  19.  Case  8  was  not  included  in  the 
combined  phase  function  because  it  is  believed  the  values  were  contaminated 
with  either  sunglint  or  shallow  water  color.  For  the  remaining  four  cases  the  red 
and  NIR  wavelengths  show  a  great  deal  of  variability  in  the  phase  function  plot. 
This  variability  could  come  from  the  fact  that  the  samples  for  the  four  cases  were 
taken  a  large  distance  from  each  other.  The  source  regions  for  the  smoke  and 
the  duration  over  water  may  account  for  a  significant  change  in  the  shape  of 
each  individual  phase  function.  This  would  account  for  a  highly  variable 
combined  phase  function.  Each  individual  phase  function  can  be  seen  in 
Appendix  A.  The  combined  phase  function  tends  to  be  decreasing  after  130°, 
while  the  model  phase  function  is  slowly  increasing.  This  illustrates  a  difference 
between  the  actual  and  assumed  aerosol  size  distribution. 
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Figure  1 9.  Phase  Function  for  Smoke  Cases  4-7,  over  the  Black  Sea,  with  NPS 

Aerosol  Model  #7  Phase  Function 
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Figure  20  shows  the  combined  phase  function  retrieved  for  three  smoke 
cases  northwest  of  the  United  Kingdom  over  the  North  Atlantic.  The  general 
trend  of  the  combined  phase  function  is  to  decrease  from  100° -120°  and  then 
nearly  flat  or  slightly  increasing  above  120°.  The  points  from  each  individual 
phase  function  fit  together  very  well  in  the  combined  phase  function  to  form  a 
cohesive  curve.  This  increases  the  confidence  level  that  the  phase  function 
displayed  is  the  true  phase  function  and  different  from  that  of  the  model.  It  would 
be  useful  to  have  more  data  points  at  higher  angles  to  see  if  the  actual  phase 
function  continues  to  closely  follow  the  shape  of  the  model  phase  function.  The 
smaller  difference  is  the  phase  function  from  that  of  the  model  indicates  that  the 
actual  aerosol  size  distribution  is  very  close  to  that  of  the  model.  Note  also  that 
model  five  was  selected  as  the  mean  model  for  this  phase  function  and  is  plotted 
on  Figure  20.  This  indicates  the  smoke  particles  have  a  size  distribution  with 
larger  particles  than  in  model  seven. 
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Figure  20.  Phase  Function  for  Smoke  Cases  9-1 1 ,  northwest  of  the  United  Kingdom, 

with  NPS  Aerosol  Model  #5  Phase  Function 
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Figure  21  shows  the  combined  phase  function  retrieved  for  three  smoke 
cases  in  a  different  area  of  the  North  Atlantic  to  the  west  of  the  United  Kingdom. 
The  general  trend  of  the  combined  phase  function  is  to  decrease  in  value  as  the 
backscatter  angle  increases.  There  appears  to  be  more  irregularities  in  the 
phase  function  than  in  the  previous  United  Kingdom  case.  The  general  shape  of 
the  phase  function  matches  fairly  well  with  the  model  phase  function. 
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Figure  21 .  Phase  Function  for  Smoke  Cases  1 2-1 4,  west  of  the  United  Kingdom,  with 

NPS  Aerosol  Model  #7  Phase  Function 


34 


Figure  22  illustrates  the  combined  phase  function  retrieved  from  four 
smoke  cases  near  mainland  Japan.  It  depicts  a  very  flat  curve  with  only  slightly 
increasing  values  at  each  end.  This  phase  function  had  the  largest  optical  depth 
values  of  all  the  smoke  cases.  It  also  had  some  of  the  highest  backscattering 
angles.  Model  number  six  indicates  that  the  smoke  particles  have  a  size 
distribution  with  larger  particles  than  in  model  seven.  The  actual  phase  function 
is  similar  in  shape  to  the  model  phase  function.  There  are  some  inconsistencies 
in  the  values  between  the  SeaWiFS  and  Terra  satellites  near  the  150°  angle. 
The  three  SeaWiFS  values  are  quite  a  bit  higher  than  the  three  Terra  values  in 
both  the  red  and  NIR  phase  functions.  The  relative  high  point  near  135°  is  also 
a  SeaWiFS  data  point.  This  also  could  be  an  indication  of  a  SeaWiFS  calibration 
issue.  Figure  23  shows  the  same  phase  function  without  the  SeaWiFS  data.  It  is 
a  much  smoother  curve  that  more  closely  matches  the  model  phase  function.  . 
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Figure  22.  Phase  Function  for  Smoke  Cases  15-18,  near  Mainland  Japan,  with  NPS 

Aerosol  Model  #6  Phase  Function 
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Japan  --  Combined  w/o  SeaWiFS 
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Figure  23.  Phase  Function  without  SeaWiFS  data,  for  Smoke  Cases  1 5-1 8,  near 
Mainland  Japan,  with  NPS  Aerosol  Model  #6  Phase  Function 


B.  DUST 

A  total  of  15  dust  cases  were  observed.  Each  case  had  four  or  five 
satellites  over  the  aerosol  within  a  two  to  four  hour  time  period.  Each  specific 
case  can  be  seen  in  detail  in  Appendix  B.  The  dust  cases  were  run  on  the  NPS 
aerosol  model  and  the  S12  values  are  displayed  in  Figure  24.  The  tighter  the 
grouping  of  the  S12  values  the  greater  the  confidence  in  the  multiple  satellite 
technique.  The  tighter  grouping  indicates  that  the  satellites  are  viewing  similar 
ratios  of  the  red/NIR  wavelengths.  All  of  the  groupings  had  small  variance 
among  the  values.  Note  that  the  NOAA-16  values  were  a  little  higher  than  the 
rest  of  the  grouping  in  cases  five  through  eight. 
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Figure  24.  S12  Ratio  for  Dust  Cases 
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NPS  Aerosol  Model  (AOD63)--Dust 
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Figure  25.  NPS  Aerosol  Model  Optical  Depth  (AOD63)  for  the  Dust  Cases 


The  aerosol  optical  depths  shown  in  Figures  25  and  26  indicate  a  large 
variability  in  the  values  when  the  optical  depths  were  large.  The  optical  depth  of 
an  aerosol  tends  to  change  more  rapidly  when  the  values  are  high.  The  aerosol 
appears  to  be  spreading  out  over  the  water  and  losing  its  constant  large  optical 
depth  value.  Another  possible  explanation  is  that  the  aerosol  particles  are 
undergoing  hydration  from  being  over  the  water  for  an  extended  period  of  time. 
In  addition,  it  was  harder  to  pinpoint  and  track  a  significant  feature  in  the  aerosol 
without  its  characteristics  changing  when  the  optical  depths  were  very  high. 
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None  of  the  sensors  seemed  to  have  a  bias  relative  to  the  other  sensors. 
The  highest  optical  depths  in  each  case  seemed  to  be  distributed  fairly  evenly 
among  the  different  sensors.  This  held  true  for  the  optical  depths  at  both  the  red 
and  NIR  wavelengths. 
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Figure  26.  NPS  Aerosol  Model  Optical  Depth  (AOD86)  for  the  Dust  Cases 
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Figure  27  reveals  a  large  variability  in  the  model  selection  for  cases  five 
through  eight.  These  are  the  dust  cases  where  the  optical  depths  were  the 
lowest  and  the  size  distribution  of  the  aerosol  lies  in  the  transition  zone  between 
models  one  and  seven.  The  model  selection  for  NOAA-16  is  the  highest  for 
these  cases  as  well  as  cases  13  and  14.  The  typical  model  selection  for  the  dust 
cases  was  model  one.  This  is  expected  since  model  number  one  has  the  largest 
particle  sizes  in  its  distribution.  As  seen  in  the  smoke  cases,  when  there  was 
disagreement  in  the  model  selections,  the  SeaWiFS  sensor  had  the  lowest  model 
value. 
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Figure  27.  NPS  Aerosol  Model  Phase  Function  Pick  for  the  1 5  Dust  Cases 
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Displayed  in  Figure  28  is  the  combined  phase  function  retrieved  for  the 
four  dust  cases  over  the  Straight  of  Hormuz  near  the  Persian  Gulf.  The  general 
trend  of  the  phase  function  is  to  be  relatively  flat  with  a  sharp  rise  after  140°. 
The  increase  in  value  after  140°  is  representative  of  what  is  seen  in  the  model 
phase  function.  There  are  some  differences  between  the  shape  of  the  actual  and 
model  phase  function.  This  indicates  that  there  are  some  differences  between 
the  actual  aerosol  size  distribution  and  the  size  distribution  in  the  model. 
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Figure  28.  Phase  Function  for  Dust  Cases  1-4,  near  the  Persian  Gulf,  with  NPS 

Aerosol  Model  #1  Phase  Function 
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Figure  29  shows  the  combined  phase  function  retrieved  for  four  dust 
cases  off  the  southern  coast  of  Pakistan.  The  combined  phase  function  plot  for 
the  red  and  NIR  wavelengths  is  very  irregular  and  does  not  follow  a  smooth 
curve.  The  optical  depths  for  these  dust  cases  were  relatively  low.  These  were 
the  cases  where  the  NPS  aerosol  model  had  a  difficult  time  also  in  determining 
an  aerosol  model  number.  The  model  number  selected  varied  greatly  with  each 
satellite  within  each  individual  case  (Figure  27).  It  is  probable  that  this 
environment  is  changing  over  the  four  hour  time  period  during  the  satellite 
passes.  With  the  low  optical  depth  values  a  small  change  in  aerosol  properties 
can  quickly  lead  to  a  relatively  large  change  in  the  model  number  selected.  The 
individual  cases  that  make  up  this  phase  function  do  not  follow  a  smooth  curve. 
Thus,  when  the  four  individual  cases  were  combined  the  result  is  an 
exceptionally  coarse  graph.  The  individual  cases  can  be  seen  in  Appendix  B.  Of 
note  is  that  model  number  three  was  the  mean  model  selected  for  these  cases. 
This  indicates  that  the  actual  particles  are  smaller  in  size  than  in  the  other  dust 
cases. 


Figure  29.  Phase  Function  for  Dust  Cases  5-8,  near  Pakistan,  with  NPS  Aerosol 

Model  #3  Phase  Function 
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Figure  30  shows  the  combined  phase  function  retrieval  for  four  dust 
cases  off  the  coast  of  Libya.  This  graph  had  some  unexpected  features.  Some 
of  the  values  seemed  extremely  high  while  others  appeared  low.  The  four  points 
that  create  a  ridge  in  the  graph  near  120°  all  came  from  the  AVHRR  sensor  on 
NOAA-17.  It  is  possible  that  these  data  values  were  contaminated  by  sunglint, 
causing  an  unusually  high  value.  Sunglint  could  not  be  ruled  out  based  on  the 
geometry  between  the  satellite  zenith,  solar  zenith,  and  their  relative  azimuth. 
The  visible  indication  of  the  sunglint  could  be  masked  by  the  heavy  dust  over  the 
water.  It  could  also  be  that  the  optical  depth  was  not  remaining  constant  from  one 
satellite  pass  to  the  next.  This  would  also  account  for  the  significant  amount  of 
fluctuation  in  the  curve  above  140°. 
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Figure  30.  Phase  Function  for  Dust  Cases  9-12,  near  Libya,  with  NPS  Aerosol  Model 

#1  Phase  Function 
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Figure  31  shows  the  combined  phase  function  retrieved  for  three  dust 
cases  over  the  Red  Sea.  The  function  is  fairly  flat  until  about  160°  then  the 
values  have  a  significant  jump  for  the  last  three  points.  From  looking  at  the 
model  phase  function  we  would  expect  a  rapid  increase  in  the  actual  phase 
function  after  140°.  There  are  some  definite  inconsistencies  in  the  shape  of  the 
actual  phase  function  compared  to  the  model  phase  function.  This  is  most  likely 
due  to  differences  in  the  aerosol  size  distribution  from  the  model  size  distribution. 
The  individual  phase  functions  for  this  case  as  well  as  the  thirteen  other  dust 
cases  can  be  seen  in  Appendix  B. 
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Figure  31 .  Phase  Function  of  Dust  Cases  1 3-15,  over  the  Red  Sea,  with  NPS 

Aerosol  Model  #1  Phase  Function 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  retrieval  of  an  aerosol  scattering  phase  function  from  multiple  polar 
orbiting  satellites  proved  to  be  a  difficult  but  attainable  task.  A  total  of  33  phase 
functions  were  derived,  18  of  which  had  smoke  aerosol  and  the  remaining  15 
were  for  dust  aerosol.  Each  of  the  cases  was  viewed  by  four  to  nine  satellite 
passes  within  a  two  to  four  hour  time  frame.  The  red  and  NIR  radiance  values 
were  used  along  with  the  aerosol  backscattering  angles  to  determine  the  phase 
function. 

The  optical  depths  necessary  for  plotting  the  phase  functions  were 
primarily  retrieved  from  the  NPS  aerosol  model.  The  optical  depths  from  the 
model  grouped  nicely  within  each  case  study.  The  greatest  variability  among  the 
satellites  occurred  when  the  optical  depths  were  the  largest.  The  Si2  ratios  also 
grouped  together  well  within  each  of  the  smoke  and  dust  cases.  The  smaller 
aerosols,  like  smoke,  had  a  larger  S12  ratio  than  the  dust  particles.  The  higher 
S12  ratio  led  to  a  higher  model  selection  for  the  smoke  cases.  The  model 
selected  most  often  was  seven  for  the  smoke  cases  and  one  for  the  dust  cases. 

The  phase  function  from  the  NPS  model  selected  was  compared  to  the 
phase  function  retrieved  from  the  multiple  satellite  technique.  Having  a  similar 
phase  function  to  that  of  the  model  indicates  the  size  distribution  used  in  the 
model  is  similar  to  that  of  the  actual  aerosol.  When  the  phase  functions  were 
different  than  that  of  the  model  it  indicated  the  aerosol  size  distribution  was 
different  than  that  of  the  model  or  that  the  aerosol  was  not  maintaining  constant 
properties  between  each  satellite  pass.  The  phase  functions  from  the  smoke 
cases  tended  to  match  the  model  phase  functions  more  closely  than  the  dust 
cases.  Although  the  phase  functions  for  the  smoke  cases  were  similar  to  those 
of  the  model,  they  were  not  identical.  This  suggests  that  the  actual  aerosol  size 
distribution  is  close  but  not  identical  to  the  model  size  distribution. 
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The  phase  functions  for  the  dust  cases  did  not  match  the  model  phase 
functions  very  well.  Some  of  the  phase  functions  appeared  very  irregular  and  did 
not  follow  a  smooth  curve;  this  implies  a  rapidly  changing  aerosol.  The  aerosol 
did  not  maintain  consistent  properties  between  each  satellite  pass.  Its  optical 
depth  was  most  likely  changing  due  to  a  spreading  out  of  the  aerosol  or  the 
hydration  of  the  aerosol  particles.  The  Pakistan  and  Libya  cases  had  the  most 
irregularities,  while  the  Persian  Gulf  and  Red  Sea  cases  had  the  least. 

Finally,  there  seems  to  be  some  inconsistencies  with  the  SeaWiFS 
sensor.  There  is  possibly  a  calibration  issue  similar  to  that  of  the  AVHRR  sensor 
on  board  the  NOAA  satellites.  Further  examination  is  required  to  conclusively 
determine  if  a  problem  exists. 


B.  RECOMMENDATIONS 

There  are  several  follow-on  recommendations  that  would  enhance  the 
value  of  this  study.  Now  that  a  representative  phase  function  of  the  actual 
aerosol  has  been  retrieved,  the  next  logical  step  would  be  to  determine  the 
aerosol  size  and  number  distribution  necessary  to  produce  that  phase  function. 
With  a  more  accurate  picture  of  the  aerosol  size  distribution  and  composition  a 
more  complete  awareness  of  the  battlespace  will  exist.  With  that  knowledge 
comes  a  better  understanding  of  what  the  impacts  will  be  to  the  warfighter. 

Another  recommendation  is  to  measure  aerosols  that  are  not  as 
concentrated  and  interrogate  areas  that  have  low  optical  depths.  Thin  dust  and 
smoke  cases  may  behave  quite  differently  than  the  thick  cases  examined  here. 
Diverse  areas  of  the  world  have  very  different  aerosol  issues.  The  yellow  wind 
from  the  Loess  Plateau  and  the  Brown  Wind  from  China  are  several  worthy  of 
study.  Determining  the  phase  function  of  the  background  environment  with  no 
aerosol  would  also  prove  useful. 
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In  the  future  the  number  of  satellites  in  orbit  will  increase  and  the  quality  of 
sensors  onboard  will  improve.  With  the  introduction  of  systems  like  the  National 
Polar-orbiting  Operational  Environmental  Satellite  System  (NPOESS)  the  number 
of  data  points  in  each  phase  function  can  be  increased.  The  more  precise  the 
phase  function  retrieval  the  more  exact  the  information  will  be  about  that  aerosol. 
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APPENDIX  A.  GRAPHS  OF  INDIVIDUAL  SMOKE  CASES 


A.  HUDSON  BAY  SMOKE  CASES 


Figure  32.  Three  smoke  cases  over  Hudson  Bay 
(from  NASA  GSFC,  2005b) 
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Figure  33.  Smoke  case  1  over  Hudson  Bay 
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Figure  34.  Smoke  case  2  over  Hudson  Bay 
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Figure  35.  Smoke  case  3  over  Hudson  Bay 


B.  BLACK  SEA  SMOKE  CASES 


Figure  36.  Four  smoke  cases  over  the  Black  Sea,  the  fifth  case  was  north  of  image 

(from  NASA  GSFC,  2005b) 
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Figure  37.  Smoke  case  4  over  the  Black  Sea 
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Figure  38.  Smoke  case  5  over  the  Black  Sea 
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Figure  39.  Smoke  case  6  over  the  Black  Sea 
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Figure  40.  Smoke  case  7  over  the  Black  Sea 
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Figure  41 .  Smoke  case  8  over  the  Black  Sea 
C.  UNITED  KINGDOM  SMOKE  CASES 


Figure  42.  Three  smoke  cases  northwest  of  the  United  Kingdom 

(from  NASA  GSFC,  2005b) 
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Figure  43.  Smoke  case  9  northwest  of  the  United  Kingdom 
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Figure  44.  Smoke  case  10  northwest  of  the  United  Kingdom 
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Figure  45.  Smoke  case  1 1  northwest  of  the  United  Kingdom 


D.  UNITED  KINGDOM  SECOND  SET  OF  SMOKE  CASES 


Thick  □  :  Med 
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Figure  46.  Three  smoke  cases  west  of  the  United  Kingdom 

(from  NASA  GSFC,  2005b) 
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Figure  47.  Smoke  case  1 2  west  of  the  United  Kingdom 
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Figure  48.  Smoke  case  1 3  west  of  the  United  Kingdom 
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Figure  49.  Smoke  case  1 4  west  of  the  United  Kingdom 
E.  JAPAN  SMOKE  CASES 

Japan 

May  23,  2003 


Figure  50.  Four  smoke  cases  near  mainland  Japan 
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Figure  51 .  Smoke  case  1 5  near  mainland  Japan 
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Figure  52.  Smoke  case  16  near  mainland  Japan 
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Figure  53.  Smoke  case  1 7  near  mainland  Japan 
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Figure  54.  Smoke  case  1 8  near  mainland  Japan 
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APPENDIX  B.  GRAPHS  OF  INDIVIDUAL  DUST  CASES 


A.  PERSIAN  GULF  DUST  CASES 


Figure  55.  Four  dust  near  the  Persian  Gulf 
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Figure  56.  Dust  case  1  near  the  Persian  Gulf 
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Figure  57.  Dust  case  2  near  the  Persian  Gulf 
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Figure  58.  Dust  case  3  near  the  Persian  Gulf 
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Figure  59.  Dust  case  4  near  the  Persian  Gulf 
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B.  PAKISTAN  DUST  CASES 


Figure  60.  Four  dust  cases  near  Pakistan 
(from  NASA  GSFC,  2005b) 
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Figure  61 .  Dust  case  5  near  Pakistan 
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Figure  62.  Dust  case  6  near  Pakistan 
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Figure  63.  Dust  case  7  near  Pakistan 
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Figure  64.  Dust  case  8  near  Pakistan 


Figure  65.  Four  dust  cases  near  Libya 
(from  NASA  GSFC,  2005b) 
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Figure  66.  Dust  case  9  near  Libya 
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Figure  67.  Dust  case  10  near  Libya 
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Figure  68.  Dust  case  1 1  near  Libya 
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Figure  69.  Dust  case  12  near  Libya 
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D.  RED  SEA  DUST  CASES 


Figure  70.  Three  dust  cases  over  the  Red  Sea 
(from  NASA  GSFC,  2005b) 
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Figure  71 .  Dust  case  1 3  over  the  Red  Sea 
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Figure  72.  Dust  case  1 4  over  the  Red  Sea 
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Figure  73.  Dust  case  1 5  over  the  Red  Sea 
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E.  NAMIBIA  DUST  CASE 

It  can  be  seen  in  Figure  74  that  there  is  dust  blowing  off  the  coast  of 
Namibia.  This  however,  is  a  case  where  the  dust  spread  out  quickly  as  it  left  the 
coast.  The  aerosol  properties  of  the  dust  plume  tracked  in  the  box  did  not  remain 
constant  for  each  satellite  image.  The  optical  depth  has  most  assuredly  changed 
as  the  dust  has  spread  out.  Therefore,  this  case  was  not  able  to  be  used  for 
evaluation  and  represents  a  specific  limitation  to  this  method. 


Figure  74.  Sequence  of  satellite  images  showing  a  dust  feature  changing  too  rapidly 

to  be  used  a  case  study 
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APPENDIX  C.  AERONET  OBSERVATIONS 


All  of  the  graphics  in  this  appendix  were  generated  at  the  NASA,  Goddard 
Space  Flight  Center  web  site  (http://aeronet.gsfc.nasa.gov/data_menu.html) 
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PI  :  Alain_Eoyer  and  Norm_0' Neill,  alain . royer&uslierbro 
Level  1.0  AOT;  Data  from  9  SEP  2003 
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Figure  75.  AERONET  data  from  Kuujjuarapik,  south  of  Fludson  Bay 
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Moldova  ,  N  47  00' 00",  E  28  48' 57" ,  Alt  205  m, 
PI  :  Brent_HoIben,  brentfraer onet.gsfc.nasa.gov 
Level  1.0  AOT;  Data  from  2G  MAR  2003 
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Figure  76.  AERONET  data  from  Moldova,  west  of  Black  Sea 


Noto  ,  N  37  20' 02",  E  137  08' 13" ,  Alt  200  m, 
PI  :  Itaru_Sano,  sano@ixifo.kindai.ac.jp 
Level  1.0  AOT;  Data  from  23  WAY  2003 
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Figure  77.  AERONET  data  from  Noto,  on  mainland  Japan 
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Dhabi  ,  N  24  28' 51",  E  54  22' 58“ Alt  15  m, 

PI  :  Brent.  Bolhen,  brent&aer onet.gsfc.nasa.gov 
Level  1.0  AOT;  Data  from  13  DEC  2003 
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Figure  78.  AERONET  data  from  Dhabi,  near  Persian  Gulf 


Lampedusa  ,  N  35  31' 01",  E  12  37' 55" ,  Alt  45  m, 

PI  :  Sergio_Pugnaghi  and  Renato_Santangelo,  pugnaghi . ser< 
Level  1.0  AOT;  Data  from  26  APR  2004 
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Figure  79.  AERONET  data  from  Lampedusa,  near  Libya,  26  Apr  2004 
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Figure  80.  AERONET  data  from  Lampedusa,  near  Libya,  27  Apr  2004 
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